Search for: All records

Abstract We present a new suite of numerical simulations of the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over $Z\prime \equiv Z/Z_{\odot }\sim 0.1-3$, gas surface density Σ_gas∼ 5–150M_⊙pc⁻², and stellar surface density Σ_star∼ 1–50M_⊙pc⁻². The range of emergent star formation rate surface density is Σ_SFR∼ 10⁻⁴–0.5M_⊙kpc⁻²yr⁻¹, and ISM total midplane pressure isP_tot/k_B= 10³–10⁶cm⁻³K, withP_totequal to the ISM weight $\mathcal{W}$. For given Σ_gasand Σ_star, we find ${\mathrm{\Sigma }}_{\mathrm{SFR}}\propto Z{\prime }^{0.3}$. We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. The total midplane pressure consists of thermal, turbulent, and magnetic stresses. We characterize feedback modulation in terms of the yield ϒ, defined as the ratio of each stress to Σ_SFR. The thermal feedback yield varies sensitively with both weight and metallicity as ${\mathrm{\Upsilon }}_{\mathrm{th}}\propto {\mathcal{W}}^{-0.46}Z{\prime }^{-0.53}$, while the combined turbulent and magnetic feedback yield shows weaker dependence ${\mathrm{\Upsilon }}_{\mathrm{turb}+\mathrm{mag}}\propto {\mathcal{W}}^{-0.22}Z{\prime }^{-0.18}$. The reduction in Σ_SFRat low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved. 

Abstract Stellar winds contain enough energy to easily disrupt the parent cloud surrounding a nascent star cluster, and for this reason they have long been considered candidates for regulating star formation. However, direct observations suggest most wind power is lost, and Lancaster et al. recently proposed that this is due to efficient mixing and cooling processes. Here we simulate star formation with wind feedback in turbulent, self-gravitating clouds, extending our previous work. Our simulations cover clouds with an initial surface density of 10 2 –10 4 M ⊙ pc −2 and show that star formation and residual gas dispersal are complete within two to eight initial cloud freefall times. The “efficiently cooled” model for stellar wind bubble evolution predicts that enough energy is lost for the bubbles to become momentum-driven; we find that this is satisfied in our simulations. We also find that wind energy losses from turbulent, radiative mixing layers dominate losses by “cloud leakage” over the timescales relevant for star formation. We show that the net star formation efficiency (SFE) in our simulations can be explained by theories that apply wind momentum to disperse cloud gas, allowing for highly inhomogeneous internal cloud structure. For very dense clouds, the SFE is similar to those observed in extreme star-forming environments. Finally, we find that, while self-pollution by wind material is insignificant in cloud conditions with moderate density (only ≲10 −4 of the stellar mass originated in winds), our simulations with conditions more typical of a super star cluster have star particles that form with as much as 1% of their mass in wind material. 

Abstract Molecular clouds are supported by turbulence and magnetic fields, but quantifying their influence on cloud life cycle and star formation efficiency (SFE) remains an open question. We perform radiation magnetohydrodynamic simulations of star-forming giant molecular clouds (GMCs) with UV radiation feedback, in which the propagation of UV radiation via ray tracing is coupled to hydrogen photochemistry. We consider 10 GMC models that vary in either initial virial parameter (1 ≤ α vir,0 ≤ 5) or dimensionless mass-to-magnetic flux ratio (0.5 ≤ μ Φ,0 ≤ 8 and ∞ ); the initial mass 10 5 M ⊙ and radius 20 pc are fixed. Each model is run with five different initial turbulence realizations. In most models, the duration of star formation and the timescale for molecular gas removal (primarily by photoevaporation) are 4–8 Myr. Both the final SFE ( ε * ) and time-averaged SFE per freefall time ( ε ff ) are reduced by strong turbulence and magnetic fields. The median ε * ranges between 2.1% and 9.5%. The median ε ff ranges between 1.0% and 8.0%, and anticorrelates with α vir,0 , in qualitative agreement with previous analytic theory and simulations. However, the time-dependent α vir ( t ) and ε ff,obs ( t ) based on instantaneous gas properties and cluster luminosity are positively correlated due to rapid evolution, making observational validation of star formation theory difficult. Our median ε ff,obs ( t ) ≈ 2% is similar to observed values. We show that the traditional virial parameter estimates the true gravitational boundedness within a factor of 2 on average, but neglect of magnetic support and velocity anisotropy can sometimes produce large departures from traditional virial parameter estimates. Magnetically subcritical GMCs are unlikely to represent sites of massive star formation given their unrealistic columnar outflows, prolonged lifetime, and low escape fraction of radiation.